Non-Abelian Yang- Mills— Higgs vortices 
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In this Letter we present new, genuinely non-Abelian vortex solutions in SU(2) Yang-Mills-Higgs 
theory with only one isovector scalar field. These non-Abelian solutions branch off their Abelian 
counterparts (Abrikosov-Nielsen-Olesen vortices) for precise values of the Higgs potential coupling 
constant A. For all values of A, their energies lie below those of the Abelian energy profiles, the 
latter being logarithmically divergent as A — s> oo. The non-Abelian branches plateau in the limit 
A — i> oo and their number increases with A, this number becoming infinite. For each vorticity, the 
gaps between the plateauing energy levels become constant. In this limit the non-Abelian vortices 
are non-interacting and are described by the self-dual vortices of the 0(3) sigma model. In the 
absence of a topological lower bound, we expect these non-Abelian vortices to be sphalerons. 

PACS numbers: 11.15.-q, ll.10.Kk, 11.15.Kc 
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Introduction 

Non-Abelian vortices on R 2 have attracted interest 
since a very long time. Nambu U pointed out that vor- 
tices of finite length in R 3 require monopoles at each 
cnd. Originally, they were proposed by Mandelstam [l| as 
flux tubes absorbed by non-Abelian ('t Hooft-Polyakov) 
monopoles at each end. In this picture the monopoles 
are bound, implying that in the dual picture where the 
duals of the monopoles are the quarks, one can describe 
confinement in QCD. 

The 't Hooft-Polyakov monopole is a topologically sta- 
ble and finite energy solution of the SU (2) Yang-Mills- 
Higgs (YMH) system on R 3 , where the Higgs field takes 
its values in the algebra, i. e., that it is an isovector, 
4> = 2 , 3 ), under 50(3) rotations. Topologically 
stable and finite energy vortex solutions of the gauged 
Higgs system on R 2 on the other hand are supported 
by the Abelian Higgs model, where the Higgs field is a 
complex scalar, ip = + i(j> 2 , i. e., it is an isovector 
(j) M = ((p 1 , (f> 2 ) under S O (2) rotations. This is the 
Abrikosov-Nielsen-Olesen (ANO) vortex 0. The field 
multiplets in the two models do not match. 

To construct a non-Abelian vortex on R 2 , it was re- 
alised by Nielsen and Olesen that it is necessary to have 
a model with more than one Higgs field. They chose @ 
two 50(3) isovector Higgs fields, each with its own sym- 
metry breaking potential and vacuum expectation value 
(VEV), but with the vacuum value of each oriented at 
different directions in isospace - in the simplest case be- 
ing orthogonal to each other. In this way the 50(3) 
gauge group is completely broken on the asymptotic cir- 
cle of R 2 , which is necessary for topological stability. 
Subsequently this construction was extended in models 
featuring distinct Higgs fields,generalising the SU (2) 
vortices of to SU (N) in These vortices, de- 

scribed as Zn vortices, are not genuinely non-Abelian 
since their flux is restricted to a single direction along 
the Cartan subalgebra. 

More recently, this problem was considered in the con- 



text of M = 2 supersymmetric QCD models by Hanay 
and Tong [13 by Auzzi et. al. [ll[ and by Eto et. al. [Hj]. 
The salient feature of these models is that they have both 
gauge and colour symmetries that are broken by the con- 
densate of the scalar fields in such a way that the unbro- 
ken subgroup results in orientational zero modes of the 
string, responsible for non-Abelian flux. 

Non-Abelian vortices have been studied intensively in 
the context of dual confinement in QCD (see the reviews 
(l3l | and [13] )• in addition to this physical application, 
they present important examples of cosmic strings [l5l . 
[l6j . relevant to cosmological phase transitions. 

In this Letter we have constructed non-Abelian vor- 
tices of a SU (2) YMH model with only one algebra val- 
ued, i. e., isovector, Higgs field. (Non-Abelian vortices 
in the Weinberg-Salam model were constructed in fl7j.) 
This model features exactly the same field multiplets, on 
R 2 , as the YMH system supporting the 't Hooft-Polyakov 
monopole on R 3 , differing from the two- Higgs models of 
and those supporting Z m vortices, and obviously from 
the SQCD models of [13, [n| . 

YMH model and non-Abelian Ansatz 

Our model on R 2 is described by the static Hamilto- 



+ (4Xy Tr -V + $ 



where $ = — ^cf> ■ 3 is the antihermitian isovector Higgs 
field, and Aj = —&Aj -a, j labeling the coordinate on 
R 2 , with a = (cr M ,cr 3 ), the Pauli matrices. The gauge 
field is defined by = d^A v — d^A v + [A^,A„] and 
the gauge covariant derivative is given by = + 
[Afj,, ■}. Note that only Aj, the magnetic components of 
the 5(7(2) connection A^ = (A ,Aj), appear in Eq. ((T|), 
since in the absence of a Chern-Simons term, the electric 
component of the connection Aq vanishes, by virtue of 
the non-Abelian Julia-Zee theorem [iH [l9|. 
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The energy of this model is not endowed with a topo- 
logical lower bound and a priori we would not expect 
the resulting vortices to be topologically stable. But 
the question of stability is more subtle than this. The 
(genuinely) non-Abelian vortices we have constructed 
numerically, present bifurcations from the correspond- 
ing Abelian profiles, on plots of their energies vs. the 
Higgs self-interaction coupling constant A. Remarkably, 
it turns out that each non-Abelian profile lies below the 
corresponding Abelian profile for all values of A and hence 
cannot be expected to decay into the Abelian vortex 
with higher energy than it has. Indeed, we show that 
in the A —> oo limit these non-Abelian vortices are de- 
seribed by the (stable) self-dual 'instantons' [13] of the 
0(3) sigma model on R 2 . However, since these Abelian 
vortices embedded in the non-Abelian theory at hand are 
known to be unstable, one cannot expect this feature of 
the non-Abelian vortices found here to imply stability. 
Furthermore, the Belavin-Polyakov 'instantons' feature 
an arbitrary scale, which indicates instability. Thus, we 
would expect that our non-Abelian vortices are in effect, 
sphalerons. The quantitativc stability anaysis will be 
carried out elsewhere. 

The radial Ansatz we use is 
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where we denote o? — cos nip a 1 + sin nip a 2 and (sx)j — 
(sin ip, — cosip) . Here {a, c, g, h} are funetions of r only 
and the integer n is the vortex number. This Ansatz, 
previously used to construct non-Abelian Chern-Simons- 
Higgs vortices 21] , is a consistent truncation of the most 
general Ansatz. 

Equations of motion 

Subject to the Ansatz Eq. |(3J), the Euler-Lagrange 
equations reduce to the following set of non-linear or- 
dinary differential equations, 



-v 2 (ah — c g) h . 



-r I — J = V (ah-cg)g, 
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(rh r ) r = -(ah - cg) a - 8v 2 A r [1 — (h 2 + g 2 )]h, 
r 
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- (ah -cg)c + 8v 2 A 2 r [1 - (h 2 + g 2 )} g , 



together with the constraint equation 
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The subseript r denotes ordinary differentiation with re- 
spect to r. The energy density now rcads 
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FIG. 1: Energy per vortex number E/n versus the Higgs 
potential coupling constant A for YMH solutions with n = 
1,2. 



the total energy E being given by E = 2tt J r£dr. 

The embedded Abelian solutions, namely the solutions 
to the embedded Abelian subsystem, correspond to the 
truncation {c = 0, g = 0}. These are the ANO vortices 
which play an important role in the classification of the 
non-Abelian vortices we have constructed. In particular 
we will study the dependence of these on the parameter 
A, so it is pertinent at this point to note that the critical 
configuration of the Abelian vortices corresponds to the 
value X = j. This is the Bogomol'nyi limit where the 
Abelian vortices do not interaet. 

Numerical results 

In order to generate vortex solutions to Eqs. ©-(d]), we 
impose boundary conditions such that the energy of the 
solutions is finite and both gauge and the Higgs field fune- 
tions are regular at the origin. The system of equations is 
solved numerically by means of a collocation mcthod for 
boundary-value ordinary differential equations, equipped 
with an adaptive mesh selection procedure. 

The only free parameters are n and A, since we fixed 
the unit of length by setting v = 1 in what follows. For 
fixed finite values of these parameters only a finite num- 
ber of regular solutions exist. There always exists one 
Abelian solution (ANO solution) for any n and A(^ 0). 
For small values of A this is the only possible solu- 
tion. However, as A increases new non-Abelian solutions 
braneh off the Abelian ones. With increasing A more and 
more non-Abelian branehes appear, their number becom- 
ing infinite for A = oo. For given n all the non-Abelian 
solutions have energy lower than that of their Abelian 
counterparts for each value of A. This braneh structure 
of the solutions is exhibited in Fig. 1 where the energy 
per vortex number, E /n, is plotted versus the constant 
A for n = 1, 2. 

In this figure we observe the first two non-Abelian 
branehes for n = 1 and the first three ones for n = 2. 
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FIG. 2: Location of the first branching points for n — 1,2, 3. 



FIG. 3: Energy per vortex number E/n versus the Higgs 
potential coupling constant A for YMH solutions with n = 3. 



The lowest non-Abelian braneh branehes off the Abelian 
solutions at A w 3.705, 0.975, for n = 1, 2, respectively. It 
is clearly seen that higher values of the vorticity n allow 
for new non-Abelian branehes at lower values of A. This 
is more explicitly shown in Fig. 2, where the locations 
of the first branching points for n = 1,2,3 are given in 
logarithmic scale. 

The structure of non-Abelian branehes may be labeled 
by the pair (n, m), n being the vorticity and the inte- 
ger m indicating the specific non-Abelian braneh for that 
vorticity n (lower m means lower energy). Notice that 
the Abelian solutions behave in a different way, which 
we will emphasize below. In faet, although Abelian so- 
lutions exist for any non-vanishing value of A, for each 
non-Abelian braneh there exists a minimal value of A, 
A(nm) (which depends on n and m), below which the 
non-Abelian braneh ceases to exist. In faet, at that min- 
imal value the non-Abelian braneh matehes the corre- 
sponding Abelian braneh for that value of n, so Aj™ 1 ^ 
corresponds to the branching points where non-Abelian 
branehes start to exist. 

A remarkable faet in Fig. 2 is that the gap between 
neighbouring branching points is roughly constant for 
each n on logarithmic scale for A. More precisely, the 
quantity n [log(l + Af„ in m+1) ) - log(l + Af„ in m) )] is roughly 
constant and independent of (n, m) . This feature be- 
comes more accurate for large A, revealing an underying 
structure in the non-Abelian seetor in the limit A — > oo. 
In faet, denoting the energy of the (n, m) non-Abelian 
solutions by -E(„, m ) = E^ nm ^{\), one observes in Fig. 1 
that for each m the energy per vortex number tends to a 
limit which does not depend on n but only on m. 

It turns out that 



hm — = 2wm 



(6) 



the energy per unit vorticity is equal to the energy of 



the unit vortex. Hence non-Abelian vortices with given 
m are non-interaeting in that limit. Fig. 3 shows this 
limit for n = 3 solutions. In that figure it is clearly seen 
that in the large A region the ratio E/(2Trn) approaches 
the integer value m that labels the non-Abelian branehes. 
One observes an infinite number of non-Abelian branehes 
for each n emerging from the logarithmically divergent 
Abelian profile, each converging to a finite limit. 

One can understand this feature as follows. We have 
verified that in this limit the contribution of the poten- 
tial term in Eq. ([1} to the energy of the non-Abelian 
vortices vanishes. (This contrasts with the correspond- 
ing situation for the vortices of the Abelian Higgs model.) 
Thus, the YMH theory supporting the non-Abelian vor- 
tices becomes a 0(3) sigma model on R 2 in this limit [26| . 
Likewise in our case the vanishing of the Higgs potential 
leads to the O (3) sigma model constraint, resulting in 
the S O (3) gauged O (3) model, which unlike in the WS 
case (22| , does not satisfy the Derrick scaling requirement 
for finite energy. To this end, we have verified that in this 
limit the contribution to the energy of the YM term TrF 2 
in Eq. ([T]) also vanishes, consistently with the Derrick 
scaling requirement, and that indeed the YM potential 
becomes a pure-gauge in this limit. Thus the only con- 
tribution comes from the TrZ?<E> 2 term in Eq. ([T]), which 
in this case reduces to Tr9<f> 2 of the scale invariant O (3) 
sigma model on R 2 . Our non-Abelian vortices in this 
limit are deseribed by the radially symmetric vorticity-ri 
subset of the non-interaeting self-dual Belavin-Polyakov 
ferromagnetic vortices [20l ]. 

The effect of non-Abelianness on YMH solutions af- 
fects not only the energy values, which become lower for 
non-Abelian solutions, but also to the way the energy 
is distributed throughout space. Both for Abelian and 
non-Abelian configurations, solutions are radial for n = 1 
(their energy density having the global maximum at the 
origin) and cireular for n > 1 (their energy density hav- 



4 



600 



400 : 



200 



Abelian: E=44.453 
non-Abelian m=3: E=42.188 
non-Abelian m=2: E=33.135 
non-Abelian m=l: E=18.120 




FIG. 4: Energy density for YMH solutions with n = 3 and 
A = 20.0. 



ing the global maximum at a finite non-vanishing value 
of r). For ring-shaped configurations (n > 1), as the 
solutions are more non-Abelian (lower values of m) the 
energy density profile spreads: the maximum is moved 
to higher values of r and its height decreases. In addi- 
tion, the value of the energy density at the origin tends to 
zero, the profile becoming more and more ring-like. This 
result is demonstrated in Fig. 4 where the energy density 



profiles of YMH solutions with n = 3 and A = 20.0 are 
shown. 

Conclusions 

As a final comment on the possible physical status of 
our solutions, we emphasise that the model on R 2 em- 
ployed here is precisely that which supports monopoles 
on R 3 . Interestingly, this YMH model on R 3 sup- 
ports also monopole-antimonopole (MA) solutions, con- 
structed in [25j . This describes a consistent picture where 
our vortices are candidates for flux tubes starting and 
ending on monopoles of opposite polarities. Our results 
are qualitatively consistent with the picture in (25|. In 
particular for vorticities n > 3, the energy density distri- 
bution in the MA configuration presents a ring shaped 
density situated on the symmetry plane (the R 2 plane 
where our vortices exist) much like the circles in Fig. 4. 

We conclude by noting that the vortices constructed 
are genuinely non-Abelian, but are not endowed with a 
topological lower bound. That leads us to expect that 
these non-Abelian solutions are unstable for finite val- 
ues of A, even though the limiting solutions, namely the 
Belavin-Polyakov vortices for A — > oo are stable. 
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